黄土山区煤层群分层开采地表位移规律分析

为了研究黄土山区地形对煤层群开采地表位移的影响,以榆林某矿为地质原型,利用数值软件模拟地表为近水平、正坡、负坡、凸面和凹面五种地形条件下煤层群开采。结果表明:黄土山区地形的变化,对煤层群开采顶板垂直应力影响不明显,对地表位移影响较大;单向坡开采,斜坡滑移引起的下沉量在坡顶与开采沉陷量形成“叠加”,在坡底形成“抵消”;在组合坡中,煤层群开采地表位移量相差很大,说明在变坡点附近地表位移较为复杂。     

郯庐断裂带肥东北部文集地区构造岩类型与构造变形分析

位于郯庐断裂带中南部的肥东地区是郯庐断裂带内韧性剪切带出露的主体区域之一。本文从肥东北部文集地区的构造岩组成及其变形入手,通过详细的野外观测、室内显微构造变形分析来确定该地区构造变形的几何学形态及其涡度、有限应变、分维度、差异应力、应变速率等参数。研究区主要出露角闪岩相肥东变质杂岩,构造岩以糜棱岩和糜棱岩化片麻岩为主,岩石变形强烈。根据极摩尔圆法得到的运动学涡度值为0.729~0.870,指示区域内的韧性剪切活动以简单剪切为主。几何学形态上,肥东北部文集地区构造整体呈一背形展出,其枢纽走向NE-SW,轴面倾向SE,物质及变形强度等方面均表现出一定的规律和对称性。研究区内动态重结晶石英颗粒边界具有统计学上的自相似性,其分维值D在1.143~1.208之间,且自背形核部向两翼,颗粒粒径增大,分维值减小。差异应力Δσ介于17.86~55.18MPa之间,应变速率ε值在1.960×10^-13~7.330×10^-12s^-1之间,且背形核部差异应力和应变速率较两翼大,呈近似对称性分布。通过对比以上构造参数可以发现,区内差异应力和应变速率表现出一定的正相关性,自核部向两翼二者均呈减小趋势;动态重结晶石英颗粒则与应变速率呈负相关性,即核部颗粒细小、两翼增大。本文从几何学、运动学以及显微构造变形分析等方面对该区构造变形展开精细化研究,借此来分析肥东北部地区郯庐断裂带的构造形态和运动学特征,这对进一步认识郯庐断裂带的演化过程及构造模型的建立有着重要的意义。     

An extended finite volume method and fixed-stress approach for modeling fluid injection–induced tensile opening in fractured reservoirs

Fluid injection–induced tensile opening is modeled using an extended finite volume method (XFVM). An embedded fracture strategy is used for the flow problem, that is, the fractures are discretized using finite volume segments without resolving the grid around them. Further, the discontinuities across fractures are modeled using special basis functions. The fracture openings due to enhanced fluid pressure and the associated shear slip due to traction free boundary condition on the fracture segments are both modeled using these special discontinuity basis functions. Mass transfer between fractures and matrix is modeled using the pressure difference. The enhancement of fracture storativity due to tensile opening leads to stronger coupling between flow and mechanics. An iterative scheme relying on the fixed-stress approach for fractures, which conserves the stress dependent terms over each iteration of the flow problem, has been introduced. Tensile opening has been simulated for single fractures embedded in two- and three-dimensional matrices. The convergence criterion for sequentially implicit fixed-stress scheme for fractures embedded in elastic media is established and has been validated numerically. Further, for 2D simulations, the effect of the matrix permeability for fracture propagation due to tensile opening has been studied.     

Reliability of mooring lines and piles for a permanently manned vessel in a tropical cyclone environment

The factors of safety for stationkeeping systems in current standards (ISO, API) are not derived or validated using reliability analyses. As the oil and gas exploration and production is breaking new boundaries, deploying new floating systems and moving into regions with harsher environments, it is of paramount importance to understand what level of reliability these new marine structures are achieving. This paper presents a reliability analysis of the mooring system of a Floating LNG (FLNG) vessel permanently moored and permanently manned offshore North West Australia in a tropical cyclone environment. The reliability analysis addresses both the mooring chain and the pile foundation. The analysis accounts for the long term characteristics of the environment, including the short term variability, in response to a given sea state and the variability and uncertainty in strength of the mooring chain and the pile. The stationkeeping system was analysed using detailed time domain simulations, capturing system non-linearities and low frequency oscillations as well as wave frequency responses and, thereby, reducing modelling uncertainties to a minimum.It is found that for the conditions modelled, neither the chain nor the pile meet a target reliability of 10⁻⁴/annum using the factors of safety commonly used in design following current ISO and API standards. New factors of safety are proposed to achieve this target reliability. For the pile design, one complicating factor is that current design standards do not explicitly define the exceedance probability that should be associated with the characteristic value of the undrained shear strength to be used in the design. It is demonstrated that the required factor of safety is crucially dependent on the definition of this characteristic value and on the level and the type of uncertainty in the soil strength profile. A recommendation is made regarding the definition of this characteristic value and the associated factor of safety. Furthermore, it is found that designing the mooring system to an environmental condition with a return period of 10,000 yr (as an Abnormal Limit State event), and setting the factor of safety to unity, meets the target reliability of 10⁻⁴/annum for the pile, if the characteristic undrained shear strength is a lower bound, defined in this paper by the 10th percentile value. For the chain however, this target reliability is not achieved.     

A finite element approach for predicting the full resistance profile of a spudcan deeply penetrating in dense sand overlying clay

This paper presents a finite element approach to calculate the full resistance profile of a spudcan deeply penetrating in dense sand overlying clay, in which a potential for an installing spudcan to experience a sudden uncontrolled punch-through failure exists. A modified Mohr-Coulomb model characterized by incorporating a four-phase variation of the mobilized strength and dilation parameters with an equivalent accumulated plastic strain is developed and tested for the overlying dense silica sand. An extended Tresca model is used for the strain softening of the underlying clay. A series of large deformation finite-element (LDFE) analyses are carried out, varying the strength and dilation parameters as well as the spudcan geometries. A fairly good performance of the present approach is verified by validating against groups of centrifuge tests data, allowing the numerical study to be extended parametrically. The four-phase variation of the mobilized strength and dilation parameters involved in the progressive failure of the upper dense sand is parametrically studied and extended to cover the range of sand relative densities that are of practical interest. Additionally, comparisons with the typical existing LDFE analyses using both simple and sophisticated constitutive models are carried out. It shows that the present approach performs fairly well to calculate the full resistance profile of a spudcan deeply penetration in both thin and thick dense sand overlying clay, especially the peak and post-peak resistance, within around 5% of the corresponding centrifuge tests results.     

基于新型高频声学多普勒流速剖面仪的河口近底部边界层观测初探

河口物质输运、能量交换与底边界层内的水动力过程密切相关,底边界层参数(如切应力、拖曳系数)的确定至关重要。挪威Nortek公司生产的新型声学多普勒流速剖面仪AD2CP相比传统ADCP具有高频、低噪的优点,可用于高频(16Hz)流速剖面观测,而被广泛应用于底边界层观测的ADV只能测量单点高频流速。本文采用AD2CP在长江口南槽最大浑浊带区域进行座底式观测,并与同步近底部三脚架上ADV的观测结果进行对比。结果表明,使用AD2CP测得的近底部平均流速与ADV的测量结果吻合良好;使用惯性耗散法计算了底切应力,基于ADV的单点高频流速数据计算结果为2.16×10~(-2)～5.69×10~(-1)N/m~2,基于AD2CP的结果为2.09×10~(-2)～4.26×10~(-1)N/m~2,二者范围大致相当。在此基础上,基于AD2CP数据计算出摩阻流速为4.55×10~(-3)～2.06×10~(-2)m/s、底拖曳系数范围为1.84×10~(-4)～2.49×10~(-3),与ADV的计算结果基本一致。此外,由于AD2CP可以获得高频的流速剖面数据,优于单点ADV,具备观测近底部边界层参数和边界层内湍流剖面的潜力。     

A stochastic method for the prediction of icebreaker bow extreme stresses

It is well established that the ship-ice interaction process is quite complex and associated ice loads on the icebreaker hull is a stochastic process. Obviously, novel accurate statistical methods and models should be developed and applied to estimate extreme bow stresses.This paper studies icebreaker bow stresses based on measured distribution of ice thickness in the Arctic Ocean on the way to and from the North Pole. Since the vessel route was carefully selected searching for easier ice conditions, the Arctic Ocean crossing was not a straight linear but a meandering path. Thus, the specific ship route data was biased with respect to general ice statistics in the region, but true with respect to the route specific ice data encountered by a ship navigating in that region. Therefore the route specific ice thickness data is directly needed for ship design and navigation analysis. It is assumed that captains are competent and knowledgeable, and therefore will select a route that provides the most favourable ice conditions.This paper contributes to study of the newest Chinese self-designed polar icebreaker, serving the purpose of enhancing icebreaker operational reliability. Finite Element Method software package ANSYS/LS-DYNA has been employed to simulate bow stress pattern for a particular icebreaker operating in the Arctic Ocean. Extreme bow stresses were estimated using Naess-Gaidai method. The latter is a first application of Naess-Gaidai method to a distribution with lower bound. Thus this paper aims at introducing an efficient method of estimating route-specific icebreaker extreme bow stresses.     

Conjugate ruptures and seismotectonic implications of the 2019 Mindanao earthquake sequence inferred from Sentinel-1 InSAR data

In 2019, four strong earthquakes of Mw>6.4 occurred successively in Mindanao, Philippines. Based on the reports from the USGS and PHIVOLCS, these earthquakes were dominated by strike-slip ruptures. Whether these earthquakes are temporally and spatially related remained unknown. We characterized the coseismic displacement fields during the earthquake sequence using an InSAR technique with Sentinel-1 SAR data. The InSAR deformation measurements convincingly reveal that the four earthquakes produced distinct coseismic displacement patterns. We estimated the source parameters of the earthquakes with a two-step inversion strategy. The optimal model suggests that the earthquake sequence resulted from the reactivation of a conjugate fault structure that involves two nearly vertical left-lateral strike-slip faults and two high-angle right-lateral strike-slip faults. We calculated Coulomb stress changes from the earthquake sequence, suggesting that the previous strong earthquakes had significant stress-encouraging effects on the following events. The regional velocities based on the GPS analysis suggest that the formation of this conjugate structure is mainly due to the westward movement of the subducting Philippine Sea Plate. This earthquake sequence provides a seismotectonic background for subsequent strong earthquakes and helps to better understand the formation mechanisms and seismotectonic implications of conjugate structure rupturing.     

Dynamic characteristics of marine soft clay under variable phase difference and initial static shear stress

Abstract

Under seismic loading, the soil layer is subjected to multidirectional cyclic shear stress with different amplitudes and frequencies because of the coupling of multiple shear waves and the soil element within a slope or behind a retaining wall is subjected to initial static shear stress before subjected to cyclic loading. Due to the complexity of seismic loading propagation, a phase difference exists between the initial static shear stress and cyclic shear stress. To investigate the influence of the phase difference and initial static shear stress on cyclic shear strain, cyclic modulus, and cyclic strength, a series of laboratory tests are performed on Wenzhou marine soft clay by multi-directional simple shear system, which can simulate the actual state better by controlling the horizontal cyclic stress in the x and y directions simultaneously. As the phase difference varies from 0° to 90°, the dynamic shear modulus increases and cyclic strain accumulation decreases with an increasing number of cycles. The shear strain increases with the initial shear stress.     

辛安煤矿1402工作面临空窄煤柱采掘响应及动态加固

随着我国煤炭资源去产能整合煤矿的增多，复采工作面临空窄煤柱采动失稳问题日益凸显，已严重制约矿井安全高效生产。为此，针对辛安煤矿复采1402工作面辅运巷道5号钻场临空窄煤柱稳定性控制的工程难题，运用数值模拟与理论分析相结合的方法，探究5号钻场临空窄煤柱稳定性采掘扰动响应特征，提出5号钻场临空窄煤柱动态注浆加固技术方案并开展现场应用和效果检验。研究结果表明:1402工作面辅运巷道掘进对5号钻场临空窄煤柱稳定性影响较小；在1402工作面回采期间，距5号钻场18~6 m范围，临空窄煤柱集中垂直应力由非对称马鞍形分布逐渐过渡为拱形分布；距5号钻场6 m时，临空窄煤柱承载叠加垂直应力超过煤体强度，塑性区完全贯通，极易破坏失稳；现场采用MP364型注浆材料及专用注浆设备对5号钻场临空窄煤柱前后5 m区域进行加固，动态注浆始终超前工作面10 m，通过深孔窥视和气体监测手段验证临空窄煤柱良好的封堵固化效果，保障了工作面安全回采，为我国整合矿井类似条件下煤柱稳定性控制提供借鉴和参考。移动阅读